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Global minima for semilinear optimal control 

problems 
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Abstract: We consider an optimal control problem subject to a semilinear el¬ 
liptic PDE together with its variational discretization. We provide a condition which 
allows to decide whether a solution of the necessary first order conditions is a global 
minimum. This condition can be explicitly evaluated at the discrete level. Further¬ 
more, we prove that if the above condition holds uniformly with respect to the dis¬ 
cretization parameter the sequence of discrete solutions converges to a global solution 
of the corresponding limit problem. Numerical examples with unique global solutions 
are presented. 
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1 Introduction 


Let us consider the following optimal control problem 

(P) min J(u ) = \\\y - yo||i, 2 (n ) + G|u||b(n) 

UkzU ad ^ " 

subject to the semilinear elliptic PDE 

-A y + <t>{y)=u in Slcl 2 , 

y = 0 on dfl, 

and the pointwise state constraints 

lla{x) < y(x) < yb(x), x <E I< c O. (1.2) 

We will formulate the precise assumptions on the data of the problem in Sec¬ 
tion 2. Since the state equation is in general nonlinear, the control problem 
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is nonconvex and there may be several solutions of the necessary first order 
conditions. These can be examined further with the help of second order con¬ 
ditions but those will only give local information and usually do not allow a 
decision on whether the given point is a global minimum of (P). It is exactly 
this question which is the starting point of our work. Assuming that we have 
an admissible control u which satisfies the necessary first order conditions we 
will formulate a condition on the adjoint variable that guarantees that ft is a 
solution of the control problem (P). This condition requires a certain L q - norm 
to be bounded by a constant that only depends on the data and that is known 
explicitly. While this approach is only of limited use at the continuous level, 
the situation is different when we apply our methods to a suitable discretisation 
(P/j) of (P). It turns out that we can obtain an analogous result for a discrete 
stationary point and the corresponding discrete adjoint state. But since now 
the discrete adjoint is available to us as a result of a numerical computation, we 
can check whether our condition is satisfied. If the answer is yes, Uh is a global 
minimum of (P^). Moreover we are able to make the connection back to the 
original control problem in that we show that a sequence of solutions of (P/J, 
that satisfy our condition uniformly in h converge to a global solution of (P). 

To the best of the author’s knowledge this is the first contribution to the 
study of uniqueness of solutions to semilinear elliptic optimal control problems. 
However, concerning the analysis, numerical treatment and implementation of 
semilinear optimal control problems many contributions can be found in the 
literature. Here we exemplarily mention the work [T] of Arada et al., [2] of 
Casas, and the work m of Neitzel et al. Further references can be found in 
[U Chapter 3], |5j, and in [3] , where the role of second order conditions in PDE 
constrained optimization is discussed. 


2 The optimal control problem (P) 


Let SI C M 2 be a bounded, convex and polygonal domain. We assume that 
<j> : R —> R. is of class C 2 and monotonically increasing. For our analysis we 
require the following structural assumption on <f>: 

Assumption 1: There exist r > 1 and M > 0 such that 

|<p(s)| < M<j)'{s)r for all sel. (2.1) 

Let us remark for later purposes that implies 

<f>'(s) < Cl (i + | s p), set, n = — T —— (2.2) 

v / r — 1 

|0(s)| < c 0 (l + |s| r °), set, r 0 = (2.3) 

Note that for a power nonlinearity of the form cj>(s) = |s| 9 ~ 2 s (q > 3) 


W'{s)\ = (q- 1 ){q - 2) | s | 9—3 = (q - 2 )(q - 1 )*=*[<//(s)]*=*, 


so that (2.11 is satisfied if we chose r = Solving this relation for q yields 
q = qrpp, which is an expression that we will encounter in our analysis below. 
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Using the theory of monotone ope rators one can show that for every u £ L 2 (Q) 
the boundary value problem ( |l.l| has a unique solution y £ Hq(Q) fl ff 2 (fl) 
which we denote by y = Q(u). Moreover, there exists a constant c > 0 such that 

l|y||i? 2 (n) < c(l + |M|i 2 (fi))- 

Next, suppose that K is a (possibly empty) compact subset of f2 and define 
Y ad := {z £ C(K) : y a (x) < z(x) < y b {x) for all x £ K}. 

Here, y a ,y b £ Co(f2) are given functions that satisfy y a (x) < y b (x),x £ K. 

We consider the semilinear optimal control problem 

(P) ™U 1 d J ( U ' ) := ^ ~ y °l^ 2 ( n ) + f H U Hi 2 (f2) 

subject to y = g{u) and y\ K £ Y ad , 

where 

U ad '■= {u £ L 2 (Q) : u a < u(x ) < u b a.e. in fl} 

and —oo < u a < u b < oo are given. By classical arguments Problem (P) has a 
solution u £ U ad - 

Remark 1 We note that the choice K = Cl also is possible, if the bounds 
satisfy the compatibility condition y a < 0 < y b on <9fl, which only requires 
minor modifications in the analysis. 


3 The variational discretization of (P) 


In this section we approximate Problem (P) using the variational discretization 
approach introduced in [7]. To this end, let 7 \ be an admissible triangulation 
of fl so that 

«= U 

TGTh 

We define the space of linear finite elements, 

Xho '■= {vh £ C(fl) : % is a linear polynomial on each T £ Th and Vh\gn = 0} 


and approximate (1.1) as follows: for a given u £ L 2 (fl), find y b . £ Xho such 
that 

/ Vi/ft, • Vv h + <j>(yh)vh dx= uv h dx \/v h £ X h0 . (3.1) 

Jfi Jfi 

Using a fixed point argument one can show that (3.11 has a unique solution 
Vh £ Xho which we denote by Gh(u). Finally, let us define 

Mh '■= {xj | Xj is a vertex of T £ Th , where T fl K ^0}. 

The variational discretization of Problem (P) now reads: 


(Ph) 


min J h (u ) := -\\y h - yo\\ 2 L i {n) + ^||u|| 2 2 (n) 
subject to y h = Q h {u), {yh{xj)) Xj eX h £ Y adi 
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where 


Y ad : = {(Zj)xjeAr h I y a (Xj) < Zj < Vb(Xj), Xj £ Ml J. 

We note that Problem (P/,) is still an infinite dimensional optimization problem 
since the controls are sought in U a d■ If a feasible point for (P^) exists, stan¬ 
dard techniques yield the existence of a solution Uh £ U a d for Problem (P^). 
The typical approach in order to find an optimum of (P/,) consists in trying to 
determine solutions of the necessary first order conditions. A formal analysis 
shows that these conditions read in our case: there exist multipliers ph £ W.o 
and Jlj el, Xj £ Mh such that 


/ VyhMv h +<j>(y h )v h dx = / u h v h dx Vv h £ X hQ , 

jq. J o 

/ Xp h • Vui + <t>'(yh)phVh dx = 

J n 

/ (■ yh~yo)vhdx+ ^2 yj v h(xj) Vv h £X h0 , 

XjeXh 

/ (p h + aUh)(u- u h ) dx > 0 \/u£Uad, 

Jn 

^ ^ dj( z j ~ Vh{ x j)) — 0 V ( z j)xj£ATh £ Y ad • 

Xj h 


(3.2) 


(3.3) 

(3.4) 

(3.5) 


Note that condition (3.41 is equivalent to the relation Uh = P[u a ,u b ] (- - )> so 


that the control variable is implicitly discretized and (3.21—(3.5 I amounts to 
solving a nonlinear finite-dimensional system. In order to state our main result 
of this section we introduce the following constant: 


r](a, r) := a^C q 


~M~ 


r — 1 \ l=r 


' 2r - 1 ; 


,i/9 r . 1 A p p/2(2_p)f-i. (3.6) 


<7 


Here, q := ir _ , p := —— while M and r appear in (2.1). Furthermore, C Q is 
‘ 1 rq 

an upper bound on the optimal constant in the Gagliardo-Nirenberg inequality 

2 q — 2 

\\f\\ L i < C\\f\\l 2 \\Xf\\ L l (q > 2). For our purposes it will be important to 
specify a constant C q that is as sharp as possible. Lemma 6.3 in the Appendix 


will give three such bounds, two of which can be found in the literature, while 
the third is new to the best of our knowledge. Let us now formulate the main 
result of this section. 


Theorem 3.1 Suppose that u h £ U a d,Vh £ X h o,Ph £ X h o, (yj) Xj eX h a 


solution of (3.2)-(3.5). If 


IMU-qn) < y(u,r), 


(3.7) 


then Uh is a global minimum for Problem (P/t). If the inequality (3.7) is strict, 
then Uh is the unique global minimum. 
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Proof: Let Uf, £ U a d be a feasible control, yh = Gh(uh ) the associated state 
with {y h {xj)) Xj eU h e Y* d . We have 

la f 

Jh{u h ) - Jhiu h ) = 2 112/ft- ^ ^111,2(0) + — \\uh — w/J^ 2 (n) + a J^u h (u h - u h ) dx 

+ (yh- yo)(yh - yh) dx (A) 

J n 


Using := y h - yh in fl3.3| ) we get 

la f 

(^) = 2 ^lli 2 (fi) + o\\ Uh - u h\\ 2 L 2 {n) + a J u h (u h -u h )dx 

+ / Vph -V(yh - Vh) + {yh)Ph(yh - yh) dx - ^ Pj(y h (xj) - y h (xj)) 

xjeAfh 

la f 

> ^\\yh-Vh\\\2 {n) +-\\ u h-u h \\L2 {Q) +aJ^u h (u h -u h )dx 

+ / ^Ph -V(yh - Vh) + {yh)Ph(yh - Vh) dx, (3.8) 

J n 


by (3.51. Using (3.11 for yh and yh with test function pi,, we get 


/ Vp/i ■ V (yh — yh) dx = / (uh — u h )ph dx — / (4>(y h ) - 4>(yh))phdx 
JQ, JQ Jo, 

= / ( U h - Uh)phdx - / ph(yh~Vh) / <t>'{tyh + (1 - t)yh) dtdx. 

J n J n Jo 


Using this in (3.8) and recalling (3.41 we arrive at 


Jh(u h ) ~ Jh(uh) (3.9) 

> \hh - VhWhw + f II u h - U h ||| 2(n) + J^(au h +Ph)(u h - u h ) dx 

- / Ph{yh-yh) / <t>'(ty h + (1 - t)y h ) - 4> (y h ) dt dx 
Jn Jo 


IQ Jo 

> 9II yn - j/h|lz,2(n) + ~^\\ u h - w/i||t2(n) — Rh{uh), 


where 


Rh(u h ) ■= / Ph(yh-yh) / (j)'{ty h + 0--t)y h ) - <j)'(y h )dtdx. 
J n Jo 


The aim is now to estimate Rh.(uh)- To begin, Lemma 6.2 implies that 
\Rh{uh)\<L r [ \ph\\yh - Vh\ 2 ( f J (ty h + (1 ~ t)y h ) dtY dx 


in 


0 


— J-j-p 


rl I 

\ph\\yh - yh\~^y J J{ty h + (1 - t)y h ) dt\y h - yh\ 2 ^j r dx, 
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, r — In (r—l)/r 

where L r = M(---) . Next, Holder’s inequality with exponents 


q = 


3r — 2 r(3r — 2) 


gr 


and 


\Rh{uh)\ ^ Lr\\Ph\\Li(Q)\\yh 11 (f2) 


<t>'(ty h + (1 - t)y h ) dt\y h - y h \ 2 dx 


r - 1 ’ 2(r - l) 2 2r - 2 
together with Lemma |6. 3 1 yield 

[f 

2r —2 4r—4 2 

< L r C q \\Ph\\Li(n)\\yh — yh\\ L 2(Qj\\V(yh — y^W^^ 

<t>'{ty h + (1 - t)y h ) dt\y h - y h | 2 dx 

n Jo 

2r — 2 2 2 

Here we also made use of the relation -(1-) = -. Applying Lemma 

r q q 

with 


0 


a ■■= [ \V(y h -y h )\ 2 dx, b := [ [ </>'(ty h +(l-t)y h ) dt\y h -y h \ 2 dx, X:=-,fi:= 

J n Jn Jo q 


we obtain 

2r — 2 4r — 4 

\Rh{ u h)\ ^ L r Cq d r \ \ph || L® (H) \ \Vh Vh\\ L 2^ 


(3.10) 


x( [ {yh - Vh )\ 2 dx + [ [ <f>'(ty h + (1 - t)y h ) dt\y h - y h \ 2 dx\ 
\J n Jq Jo J 


where 


d r = g -i/ 9 r -i/r / 9 -P j p = 


r + q 
rq 


Using again (3.11 for yh,yh , this time with test function — yh yields 

[ NiVh -y h )\ 2 dx+ [ [ ft (ty h + {1 - t)y h ) dt\y h - y h \ 2 dx 

J £1 J £1 J 0 

< ||Uh - u h \\ L 2 (n)\\yh - y h \\mn)- 


Inserting this estimate into (3.101 and observing that + p = 2 — p we deduce 


\Rh{Uh)\ < L r Cq ’ dr||Pfc||z,®(fi)||y/i yh\\ L 2(n)W Uh Uh \\ P L 2 (£l) 

'1, 


2a 2 L r C q d r \\ph\\Li(Q.) (2 Wdh ^lli 2 (n)) ( K 2 ^ Ufl Uft lli 2 (n)) 


Applying again Lemma |6.1| this time with the choices 
1,, ,,o . ol 


a T^WVh - yhWl^n), 0 := -|| u h - u h \\~ L 2 {n) , X := 1 - n := 


_ v ._, 6 := ttII^ “ W 2 

we obtain 

_ £ 2r-2 

^ 2o; 2 L r Cq d r e r \\ph ||l9(q) 


( 9 II Vh — ?/?i|li 2 (n) + ~^W Uh — M /illi 2 (n)) > 


(3.11) 
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where 


= (i- p 2 ) 1 ~H p 2 ) § - 


> 


.)( 


1 - 2a 2 L r C a 


Using ( |3.11[ ) in ( |3.9[ ) we get 

dh {.'U-h ) Jhi^h) 

2 \\yh ~ y h \\ 2 L 2 (n) + 2 \\ Uh ~~ Uh \\~L 2 (^)j 
so that Jh{uh ) > Jh(uh) provided that 

2r — 2 

By direct calculations, we have 

2 d r e r = q~ 1/q r- 1/r p- p/2 { 2 - p) 1 " 2 . 

Hence, using the above result and the value of L r from Lemma |6.2| we can 
rewrite (3.121 as 


L 9 (fi) ~ (2o< 2 L r C q r d r e r ^ 


-l 


d r e r \\ph \\li (n)j 
(3.12) 


\\Ph\\Ll(Q) < a 2 Cq 
which is the desired result. 


1 M- 1 ( ^ g 17 V V /2 ( 2 - p) ^ 


Since the adjoint state ph and the quantity r](a, r) can be computed explicitly , 
Theorem |3.1| allows us to decide whether a function Uh which satisfies the nec¬ 
essary conditions of first order is a global minimum of (P/,). A natural question 
then is, whether a sequence {uh)o<h<h 0 of minima satisfying (3.7) uniformly in 
h converges to a global minimum of (P). We shall address this problem in the 
next section and it will be useful to have a continuous analogue of Theorem |3.1| 
A function u £ U a d satisfies the necessary first order conditions for problem (P) 
if there exist p £ L 2 (fl) and a measure p £ A4(K) such that 


/ Vy • Vv + 4>{y)v dx 

= 

f uvdx Wv£Hq(H), 

(3.13) 

J n 



Jn 


/ p(—Av) + (jf {y)pv dx 

In 

= 



[ (y- 

S> 

o 

S* 

1 

+ 

[ vdp Vu £ Hq (0) n H 2 (H), 

(3.14) 

Jn 



Ik 


/ {p+au){u 
Jn 

— u) dx 

> 

0 \/ U £ U a d, 

(3.15) 


-y)dp 

< 

0 Vz£Y ad . 

(3.16) 


IK 


It is well-known that the function p then belongs to 1 Uq’ s (H) for any 1 < s < 2 
and hence to L q [fH) for any q < oo (recall that H C M 2 ). Arguing in almost the 
same way as in the proof of Theorem |3.1| we obtain: 


Theorem 3.2 Suppose that u £ U a d,y £ Ho(fl),p e T 2 (H),p £ M(K) is 
solution of {3.13)-(3.16). If 


lblU«(n) < V(a,r), 


(3.17) 


then u is a global minimum for Problem (P). If the inequality (3.11) is strict, 
then u is the unique global minimum. 
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4 Convergence analysis 


Let (Th)o<h<h 0 be a quasiuniform sequence of triangulations of Cl. We consider 
the corresponding sequence of control problems (P^) and suppose that Uh £ U a d 
satisfies the hypotheses of Theorem |3.1| uniformly in 0 < h < ho- Thus there 


(4.1) 


exist ph £ Xho and (jij) Xj eM h satisfying (3.2)-(3.5) as well as 


U«(fi) < V(a,r), 0 <h<h 0 . 

It is convenient to introduce the measure jlh £ A4 (S4) by 


hh 'y ' ftjfixji 

Xj GAfh 


where 5 Xj is the Dirac measure at Xj. Since K C Cl, dist(a; :) -, K) < h , Xj £ 7V) ( and 
y a {x) < yb{x),x £ K there exists a compact set K c Cl, 6 > 0 and 0 < hi < h 0 
such that I\ C K and 


A f h C K, 0<h<h lt (4.2) 

Va{x) + 5 < lj(y a (x) + Vb(x)) < yb{x) - <5, x £ K. (4.3) 

Applying a smoothing procedure to x - ( y a +Vb ) £ Co(fi) we obtain a function 

w £ C™(Cl) such that 

. . S . . . . 6 

y a {x) + - < w{x) < y b {x) - x £ K. 

Let us denote by Rh : LAg(D) —> X^o the Ritz projection defined by 


XRhW-Xvhdx= / Xw-Xvhdx £ Xho- 


(4.4) 


Since RhW —¥ w uniformly in Cl, we may assume after choosing hi smaller if 
necessary that 


, . 6 . . . . 5 

y a {x) + - < R h w{x) < y b {x) - 


x £ K. 


(4.5) 


Our first step in the convergence analysis are uniform bounds on the optimal 
control Uh as well as on y b = Gh(uh) and jlh- 


Lemma 4.1 Let u h £ U ad ,y h ,Ph £ X h0 and {nj) Xj ^N h be a solution of \3.ty - 
(3.5). Then there exists a constant C > 0, which is independent of h, such 
that 

||%IU 2 (n)> ||y7i||jri(n), \\hh\\ M (k) < c - 


Proof: To begin, fix a function uq £ U a d- Inserting uq into (3.4) we infer 

u 0 (au h +ph)dx - / UhPhdx 
! Jft 


a \\ u h\\L 2 (n) < 








< lbolb 2 (n)(a|b>tlU 2 (fi) + |bbb 2 (fi)) + lb^lb 2 (Q)lb?ilb 2 (n)- 


3r - 2 

Since q = - > 3 we deduce with the help of (4.11 

r — 1 


Testing 

infer 


Uh\\L 2 (n) < C(lbo||L 2 (n) + \\Ph\\L 2 (n)) < C(lbol|L 2 (n) + lb/»IU«(fi)) < C. 

with yh, using the monotonicity of cf> and Poincare’s inequality we 


h 1 (n) < C(l + lb/j||i,2(n)) < C. 


(4.6) 


Furthermore, (2.21, (2.31 along with the continuous embedding if 1 (S2) £*(0) 

for all 1 < t < oo yield 


H0(2/fe)IU 2 (O)! lb , (y/t)IU 2 (n) < C. 


(4.7) 


In order to verify the uniform boundedness of IbhHjvqx) we first observe that 
(3.5) implies 


y h{x . )= {Vb(*j\’ ^> 0 , 

V 3 \ Va(Xj), if Hj < 0. 


As a result we deduce with the help of ( |4.5| ) 

S S 

l\\M\ M (k) = J E l^'l- E PASh(xj) - Rhw(xj)). 


xj&Nh 


Xj 


Using Vh = yh — RhW in (3.31 we may continue 


^ll Ph\\M{K) — 


Vpft, ■ '\7yh.dx — / \7p h -X7R h wdx 
1 Jn 

/ <t>'{yh)Ph{yh ~ Rhw)dx - / (y h - y 0 )(y h - R h w)dx 

Jn Jn 


= E 5 - 


(4.8) 


If we let Vh = Ph in (3.21 we obtain with the help of (14.71) and (14.1 


\Si\ = \ (u h - (j){y h ))phdx | < (Hu^Hl 2 ^) + ||<Kyblb 2 (n))IMU 2 (fi) < C. 

Jn 

Next, the definition of the Ritz projection and integration by parts yields 
S 2 = — Vph • S7wdx = / phAwdx 

Jn Jn 

so that 

I£ 2 1 < |bi»lb 2 (n)||Atu|b 2 (n) < C. 

Holder’s inequality along with (14.71), (|4.1|) and (|4.6| implies that 


bs| < W{Vh)\\L*m\\Ph\\L'<n)\\Vh - < C\\Vh - Rhw\\m(n) < C. 

Finally, 

I£ 4 1 < (lb/t|U 2 (n) + lbo||z, 2 (fi)) (lbft||z, 2 (n) + ||-R/iHU 2 (n)) < C. 
Inserting the above estimates into (4.8 1 yields the bound on |bbl.A/i(if)- ® 


9 










Now, we are in position to formulate the main theorem in this section: 


Theorem 4.2 Suppose that (v,h, UhiPhi A t h)o<h<ft 1 Is a sequence satisfying (3.2)- 
(3.5) as well as (f.l). Then 

Uh —>• u in L 2 (£l) for a subsequence h —> 0, 

where u is a global minimum for Problem (P). If 


IKIU^a) < Kr](a,r), 0 < h < hi, (4.9) 

for some 0 < n < 1, then u is the unique global solution of (P) and the whole 
sequence (uh)o<h<hi converges to u. 


Proof: From Lemma [~i~T| we deduce the existence of a subsequence h —»• 0 and 


u G L 2 (fl), y € Hq(£ 1), p G p G JA(K) such that 

u h —^ u in L 2 (fl), (4-10) 

Vh V in Hq(Q) and p h -t y in £‘(12), 1 < t < oo. (4.11) 
p h ^ P in M{K), (4.12) 

Ph -*• P in L q (n), (4-13) 


Our aim is to show that ( u, y,p,p ) is a solution of pdTsj ) -( |3.16| ). It is easy to 
see that u G U a d and that p = G(u), so that ( 3.13| is satisfied. Furthermore, 
the fact that dist(x j,K) < h,Xj G TV),, implies that supp (p) C K. Combining 
this with the bound ||M/i|lx(ff) — C we infer that 


z dp h 


zdp as h —> 0 


(4-14) 


< K 


' K 


for every sequence ( z h ) 0< h<h 1 C C(K) converging uniformly to 2 on K. Next, 
we claim that 

ph —* y uniformly in Q. (4-15) 

To see this, denote by y h G H 2 ( 17) (~l H { \ (O) the solution of 

-Ay h = u h — 4>{ph) in f2, j/ 1 = 0 on <912. 


We deduce from Lemma 


4.1 


and (4.7) that ( y fl )o<h<h 1 is bounded in H Z (Q), so 


that there exists a further subsequence and a function y G A 2 (12) ni4g(12) with 


y in H 2 (H), y h —► y in C(12). 


Since Uh — 4>{ph) — k u, — </>(y) in L 2 (12) we find that —A y = —Ay a.e. in 12. 
Hence y = p and y h —> y in (7(12). On the other hand, the definition of y h 
implies that ph = Rhy h , so that standard interpolation and inverse estimates 
imply 


II yh - j/IU°°(n) 


< \\Rhy h - y h \\L°°(n) + || y h - y||L~(f2) 

< Ch\\y h \\ H 2 (Q) + ||y ,! — 2 /||loo(q) —> 0 as h —> 0, 


since \\y h \\H 2 (n) < C. This proves (4.15l. 

Let us check that p\x G Y ac [. For a fixed point x G K we can choose a sequence 
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i x jh)o<h<hi such that x jh £ N h and \x jh - x\ < h. Since y a (xj h ) < yh{x jh ) < 
Vb(xj h ) we obtain y a (x) < y(x) < yb{x) by passing to the limit h —> 0 and using 

I fhTsl i. 

Next, let us fix z £ Y a( i and extend 2 : t o a f u nction 5 £ Q(K) satisfying y a {x) < 
z(x) < yb{x),x £ K. We obtain from (3.51, (4.141 and (4.151 

0 > Y] yj(z{xj) — y h (xj)) = / (5 — y h )dy h —>■ / (z - y)dp, 

_ r— \r Jit JK 


i&Mh 


which yields (3.16l. 

In order to derive (3.14) we fix v £ H 2 (fl) n Hq(Q) and insert Vh = RhV into 
( [373] ), i.e. 

/ \7p h ■ VRhV + 4> [xjh)PhR h v dx = / (y h - y 0 )R h vdx + / R h v dp h . 

J q J n Jit 

Using the definition of Rh and integration by parts we may write 

/ X7ph ■ VRhV dx = / V ph ■ V v dx = / ph(—Av)dx 

Jq Jq Jq 


so that (3.14) follows from passing to the limit h —► 0 taking into account (4.131, 
( f~U5j ) and Q4.14| ). 

Our next goal is to show that tfy —» ?2 in L 2 (S1). Inserting u into (3.41 and 
rearranging we infer 

“ll^lll^n) ^ / u(auh + ph)dx — / u h p h dx. (4-16) 

Jq. Jq 

The second integral can be rewritten with the help of (3.21 and (3.3), namely 
/ u h p h dx = / Vi/ft ■ 'S/phdx + / 4>(y h )phdx 

= (yh)PhVhdx + / {ph-yo)yhdx+ / y h dp h + / <j)(y h )Phdx. 

Jq Jq Jk Jq 

This relation allows us to pass to the limit in a similar way as above to give 

/ Uhphdx -> - / <f>\y)pydx + (y - yo)ydx + / ydp + / fi{y)pdx 
Jq Jq Jq Jk Jq 

= / (—A y)pdx + / <p(y)pdx = / updx , 

J Q ■/ Q n 


where we used (3.141 and the fact that y = Q(u). We can now pass to the limit 
in (4.161 and deduce that 

limsup||%||| 2(n) < ||u||i 2 (n) . 

h—> 0 

Since ||u||| 2 (n) < liindnf ||fi ft ||| 2(n) we infer that \\u h \\ L 2 (U) ||w|| L 2 (n) , which 

together with the fact that Uh —*■ u in L 2 (fl) implies that tfy —> u in L 2 (fl). 
Combining this with the weak convergence ph —*• p in L 2 (fl), one can pass to 
the limit in (3.4) to obtain 


/ (p + au)(u — u) dx > 0 \/u £ U a d, 
Jq 


(4.17) 
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which is (3.151. In conclusion we see that (u, y,p, fi) is a solution of (3.13 (-(3.16 1 . 


Furthermore, the lower semicontinuity of the L q -norm implies that 
||p||z,«(fi) ^ ii^mf \\Ph\\L«(n) < v(a,r) 


and we infer from Theorem |3.2| that u is a global minimum of Problem (P). If 
(4.9) holds, then p satisfies ||p||i 9 (n) < Krj{a,r ) < rj(a,r) and u is the unique 
global minimum of (P). A standard argument then shows that the whole se¬ 
quence (■u/ l )o</i<fc,i converges to u. ■ 


Before we go to the numerical examples, we make the following general remarks. 

Remark 2 

1. We do not require a constraint qualification such as a Slater condition to 


deduce that (u,y,p,/j,) satisfies the system (3.13)-(3.16), which represents the 
first order necessary optimality conditions for Problem (P). 


2. It is well known that ( |3.2[ )-( |3.5| can be rewritten equivalently as a system 
of semi-smooth equations and thus can be solved by a semi-smooth Newton 
method, see for instance 0 , ®, m- In particular, we can avoid the use of 
relaxation methods such as Moreau-Yosida relaxation, interior point methods, 
or Lavrentiev-type regularization. 


3. Since we solve (3.2)-(3.5) in practice on the computer, we consider Uh a global 


minimum if the inequality (3.7) is satisfied up to machine precision. Here, the 
integral ||ph||.L<i on the left hand side of this inequality is assumed to be calcu¬ 
lated exactly. However, this assumption can be achieved easily whenever q is 
an integer because in this case the function \ph\ q restricted to every triangle in 
the mesh is a (possibly piecewise) polynomial of degree q. Hence, one can use 
an appropriate quadrature rule to evaluate such an integral exactly. 


5 Numerical Examples 


In this section we consider variational discretization of the optimal control prob¬ 
lem (P) for different choices of the nonlinearity <f> and the data yo^u a , Ub , 2/cn 2/f>> cc, 
while := (0,1) x (0,1) is kept fixed in all considered examples. For the desired 
state ?/o we consider the following two choices 

A1 : 2 / 0 ( 2 ) := 2 sin(27ra;i) sin(27ra2)j 

A2 : 2 / 0 ( 2 ) := 60 + 160 ( 2 : 1 ( 2:1 — 1) + 2 : 2(22 — 1)). 


We note that in choice A 1 the desired state 2/0 vanishes on the boundary dil of 
the domain, while in choice A2 it doesn’t, see Figure [T] The numerical solution 
of the corresponding systems ( |3.2[ )-( 3.5 I is performed with the semismooth New¬ 
ton method proposed in [4j, whose extension to the treatment of finite element 
approximations of semilinear PDEs ist straightforward. All the computations 
are performed on a uniform triangulation of Cl with mesh size h = 2 —5 -\/ 2 . 


12 























Example 1 In this example we define (f>(s) := s 3 . It is easy to see that this 
nonlinearity satisfies Assumption 1 with r = 2 and M = 2x/3. Hence, in view of 


Theorem 3.1 we have q = 4 and a control Uh obtained from solving (3.2 (-(3.5) 


is a global minimum if the associated adjoint state ph satisfies 


II Ph 


||z,4(n) < 5 *2%y/2C i 'a*, 


where C, 


-l 


1.543145399297809 is the constant from Lemma 6.3 For this 


example we consider the following three cases. Let us abbreviate 
rj(a ) := r](a,2) = 5 — ®3§ V / 2C , j" 1 a5. 


Case 1 (unconstrained problem) In this case we set 


Ub — U a — OO, 

Vb = -y a = oo. 


In Table [l] we provide the values of ||ph||£ 4 , 77 (a) and J(uh) for different values 
of a where we consider the choice A1 for the desired state z/o ■ The findings are 
illustrated graphically in Figure [2] We see that for all values of a we can claim 
that v,h is a global minimum since ||ph||z,4 is less than its corresponding 77 (a). 
On the other hand, if we consider the choice A2 for 7/0 we can claim Uh is a 
global minimum only for approximately a greater than 10 -2 as it can be seen 
from Figure [3] The numerical values are provided in Table [2] 

Case 2 (constrained control) In this case we consider constraints only on the 
control, we set 


u a = - 5 , 

Ub = 5, 

Vb = -Va = 00 . 

Table [ 3 ] shows the values of ||p/t||i 4 , 77 (a) and J(uh) computed for different 
values of a while considering the choice A1 for y 0 . The graphical illustration of 
these findings are shown in Figure [4] We see that Uh is a global minimum for 
a approximately greater than 1CP 5 . The numerical results associated with the 
choice A2 are given in Table [4] and illustrated in Figure [5] In this case Uh is a 
global minimum for a approximately greater than 10 _1 . 

Case 3 (constrained state) In this case we consider constrains only on the state, 
we set 


Ub — U a — OO, 

Va = “I, 

Vb = 1- 

The numerical findings associated with choice A1 are provided in Table [5] and 
illustrated in Figure l6] For the choice A2 the results are given in Table [6] and 
illustrated in Figure [7] In both cases we see that Uh is a global minimum for all 
available values of a. 
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vq yo 




(a) yo choice Al. (b) yo choice A2. 

Figure 1 The desired state yo choices Al and A2. 


Table 1 Example [I] Case 1 with choice Al for yo'. The values of ||ph||i4, 77(a) and 
J{uh) for different values of a. 


a 

IKIU 4 

?y(a) 

J i^h) 

1.0e-06 

9.990654861172e-05 

6.776197632762e-03 

3.344560044987e-03 

1.0e-05 

9.328604940252e-04 

1.606889689070e-02 

3.128947575776e-02 

1.0e-04 

5.916313713912e-03 

3.810535956559e-02 

1.967337721757e-01 

1.0e-03 

1.322797500856e-02 

9.036204771862e-02 

4.320833160546e-01 

1 . 0 e -02 

1.509224717529e-02 

2.142821839497e-01 

4.922544738762e-01 

1 . 0 e -01 

1.530600543072e-02 

5.081431366100e-01 

4.992144829702e-01 

1 . 0 e +00 

1.532768796263e-02 

1.204997272869e+00 

4.999213370332e-01 

1 . 0 e +01 

1.532985932323e-02 

2.857498848277e+00 

4.999921325890e-01 

1 . 0 e +02 

1.533007649041e-02 

6.776197632762e+00 

4.999992132478e-01 

1.0e+03 

1.533009820744e-02 

1.606889689070e+01 

4.999999213247e-01 


Table 2 Example [I] Case 1 with choice A2 for yo'. The values of ||ph||i 4 , 77(a) and 
J{u h ) for different values of a. 


a 

\\Phh* 

77 (a) 

J(u h ) 

1.0e-06 

7.823778739727e-03 

6.776197632762e-03 

7.227759688190e+01 

1.0e-05 

2.234541300612e-02 

1.606889689070e-02 

1.065710637346e+02 

1.0e-04 

5.805844706415e-02 

3.810535956559e-02 

1.386316936362e+02 

1.0e-03 

1.125576598202e-01 

9.036204771862e-02 

1.568821491955e+02 

1 . 0 e -02 

2.290137136719e-01 

2.142821839497e-01 

1.625724420922e+02 

1 . 0 e -01 

2.997603240217e-01 

5.081431366100e-01 

1.642031427088e+02 

1 . 0 e +00 

3.061090377257e-01 

1.204997272869e+00 

1.644198126030e+02 

1 . 0 e +01 

3.066635772733e-01 

2.857498848277e+00 

1.644419766418e+02 

1 . 0 e +02 

3.067181181971e-01 

6.776197632762e+00 

1.644441976184e+02 

1.0e+03 

3.067235630566e-01 

1.606889689070e+01 

1.644444197614e+02 
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© - r](a) 


r 



10 10 


10“ 6 10“ 4 10' 2 10° 10 2 10 4 


( a ) WphWl* and »?(«) vs. «■ 


(b) J(uh) vs. a. 



(c) The optimal state yh- 


(d) The adjoint state p^. 



(e) The optimal control Uh- 

Figure 2 Example [l] Case 1 with choice A1 for yo: The values of Uphill, rj{a) and 
J(u h ) vs. a. The optimal state yn, the optimal control Uh and the adjoint state ph 
for a = 10 _1 . 
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10 10 



( a ) \\Ph\\l* and v( a ) vs. «■ 


(b) J(uh) vs. a. 



(c) The optimal state yh- 


(d) The adjoint state p^- 



(e) The optimal control Uh- 

Figure 3 Example [l] Case 1 with choice A2 for yo: The values of ||ph|| £ 4, 77(a) and 
J{uh) vs. a. The optimal state yn, the optimal control Uh and the adjoint state ph 
for a = 10 _1 . 
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Table 3 Example [l] Case 2 with choice A1 for yo: The values of ||ph||i4, 77(a) an d 
J{u h ) for different values of a. 


a 

\\Ph\W 

77 (a) 

J i^h) 

1.0e-06 

1.455724773650e-02 

6.776197632762e-03 

4.507886038196e-01 

1.0e-05 

1.455724403855e-02 

1.606889689070e-02 

4.508916148391e-01 

1.0e-04 

1.455717724977e-02 

3.810535956559e-02 

4.519082323790e-01 

1.0e-03 

1.457338622672e-02 

9.036204771862e-02 

4.612690393001e-01 

1.0e-02 

1.509224717529e-02 

2.142821839497e-01 

4.922544738762e-01 

1.0e-01 

1.530600543072e-02 

5.081431366100e-01 

4.992144829702e-01 

1.0e+00 

1.532768796263e-02 

1.204997272869e+00 

4.999213370332e-01 

1.0e+01 

1.532985932323e-02 

2.857498848277e+00 

4.999921325890e-01 

1.0e+02 

1.533007649041e-02 

6.776197632762e+00 

4.999992132478e-01 

1.0e+03 

1.533009820744e-02 

1.606889689070e+01 

4.999999213247e-01 


Table 4 Example [I] Case 2 with choice A2 for yo: The values of ||ph||i4, 77(a) and 
J(u h ) for different values of a. 


a 

\\Phh* 

?y(a) 

J{u h ) 

1.0e-06 

2.954513493743e-01 

6.776197632762e-03 

1.636849171437e+02 

1.0e-05 

2.954513728927e-01 

1.606889689070e-02 

1.636850204333e+02 

1.0e-04 

2.954526968135e-01 

3.810535956559e-02 

1.636860509190e+02 

1.0e-03 

2.954464960067e-01 

9.036204771862e-02 

1.636961799251e+02 

1.0e-02 

2.955530339094e-01 

2.142821839497e-01 

1.637871978058e+02 

1.0e-01 

2.998739300063e-01 

5.081431366100e-01 

1.642034478360e+02 

1.0e+00 

3.061090377257e-01 

1.204997272869e+00 

1.644198126030e+02 

1.0e+01 

3.066635772733e-01 

2.857498848277e+00 

1.644419766418e+02 

1.0e+02 

3.067181181971e-01 

6.776197632762e+00 

1.644441976184e+02 

1.0e+03 

3.067235630566e-01 

1.606889689070e+01 

1.644444197614e+02 


Table 5 Example [I] Case 3 with choice A1 for yo: The values of ||p/i||x,4, 77(a) and 
J{u h ) for different values of a. 


a 

\\Ph\W 

77 (a) 

J{u h ) 

1.0e-06 

1.166321621310e-04 

6.776197632762e-03 

6.248613075636e-02 

1.0e-05 

8.045399583166e-04 

1.606889689070e-02 

8.942494600427e-02 

1.0e-04 

5.009426247692e-03 

3.810535956559e-02 

2.037409649052e-01 

1.0e-03 

1.322797500856e-02 

9.036204771862e-02 

4.320833160546e-01 

1 . 0 e -02 

1.509224717529e-02 

2.142821839497e-01 

4.922544738762e-01 

1 . 0 e -01 

1.530600543072e-02 

5.081431366100e-01 

4.992144829702e-01 

1 . 0 e +00 

1.532768796263e-02 

1.204997272869e+00 

4.999213370332e-01 

1 . 0 e +01 

1.532985932323e-02 

2.857498848277e+00 

4.999921325890e-01 

1 . 0 e +02 

1.533007649041e-02 

6.776197632762e+00 

4.999992132478e-01 

1.0e+03 

1.533009820744e-02 

1.606889689070e+01 

4.999999213247e-01 
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-J{uh) 


(a) WphWl^ and v( a ) vs. «■ 


(b) J(uh) vs. a. 



(c) The optimal state yh- 


(d) The adjoint state p^. 



(e) The optimal control Uh- 

Figure 4 Example [l] Case 2 with choice A1 for yo: The values of ||ph|| £ 4 , 77 (a) and 
J{uh) vs. a. The optimal state yn, the optimal control Uh and the adjoint state ph 
for a = 10 _1 . 
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(b) J(u h ) vs. a. 



(c) The optimal state yh- 
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(e) The optimal control Uh- 


(f) The control active set (Uh = —5 inside 
the polygonal region). 


Figure 5 Example [l] Case 2 with choice A2 for yo: The values of ||ph ||_L 4 , 77 (a) and 
J{u h ) vs. a. The optimal state yh, the optimal control Uh, the control active set and 
the adjoint state ph for a = 10 -1 
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(c) The optimal state yh- 


(d) The adjoint state p^. 



(e) The optimal control Uh- 

Figure 6 Example [l] Case 3 with choice A1 for yo: The values of ||ph|| £ 4 , 77 (a) and 
J{uh) vs. a. The optimal state yn, the optimal control Uh and the adjoint state ph 
for a = 10 _1 . 
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( a ) \\Ph\\l* and v( a ) vs. «■ 




- *- J(u h ) 


(b) J(uh) vs. a. 



(c) The optimal state yh- 


(d) The adjoint state p^- 



(e) The optimal control Uh- 

Figure 7 Example [l] Case 3 with choice A2 for yo: The values of ||ph|| £ 4, 77(a) and 
J{uh) vs. a. The optimal state yn, the optimal control Uh and the adjoint state ph 
for a = 10 _1 . 
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Table 6 Example [l] Case 3 with choice A2 for yo : The values of ||z,4, 77 (a) and 
J(u h ) for different values of a. 


a 

\\Ph\W 

77 (a) 

J(Uh) 

1.0e-06 

8.727496956489e-04 

6.776197632762e-03 

1.525635329141e+02 

1.0e-05 

6.303449080470e-03 

1.606889689070e-02 

1.536018906075e+02 

1.0e-04 

2.143214405409e-02 

3.810535956559e-02 

1.559131574621e+02 

1.0e-03 

8.541044896637e-02 

9.036204771862e-02 

1.600259053817e+02 

1.0e-02 

1.596641521237e-01 

2.142821839497e-01 

1.627648901943e+02 

1.0e-01 

2.997603240217e-01 

5.081431366100e-01 

1.642031427088e+02 

1.0e+00 

3.061090377257e-01 

1.204997272869e+00 

1.644198126030e+02 

1.0e+01 

3.066635772733e-01 

2.857498848277e+00 

1.644419766418e+02 

1.0e+02 

3.067181181971e-01 

6.776197632762e+00 

1.644441976184e+02 

1.0e+03 

3.067235630566e-01 

1.606889689070e+01 

1.644444197614e+02 


Case 4 The following example is taken from cn Section 7]. In particular, 
= s 3 and 

Ub — U a — OO, 

Vb = oo, 

, , 2 {1, , 1, 
y a {x) = -- +mm I -{X! +X 2 ),-{l + x 1 

Vo = -1 
a = 10” 3 . 

The numerical findings for this case are given in Table [7] and they are 
illustrated graphically in Figure [8] It is clear that Uh is a global minimum for 
the given values of a. 


— X2), ^(1 — *1 + X2), 1 — i(a;i+x 2 ) ), 


22 





9 ^ 



0.5 

-T-T--*-♦-♦- 

JD 


/ 

XT ' 

0.45 

* 

JD' 


. 


0.4 

/ 

_ 0" 


/ 

„ 0 " 


4 

^0' 

IS 0.35 

/ 


h 's 

' 






0.3 


*' - e - 17 (a) 





0.25 

/ 



& 



e 


.4 10 ' 6 10 “ 4 10 2 10° 10 2 10 4 
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(e) The optimal control Uh- (f) The multiplier 

Figure 8 Example [I] Case 4: The values of ||ph||z,4, 7j(a) and J(uh) vs. a. The 
optimal state yn, the optimal control Uh, the adjoint state ph and the multiplier 
for a = 10~ 3 . 
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Example 2 In this example we define (f>(s) := s 5 . We see that Assumption 
1 is satisfied with 

4 , ,, 20 

r = - and M =» . 

3 5* 

Hence, in view of Theorem |3.1| we have q = 6 and a control Uh obtained from 
solving (3.2)-(3.5) is a global minimum if the associated adjoint state ph satisfies 

1124 


\\Ph\\L»(n) < 


13h2bV3 C ' 6 


where C fi 


1.271251384316953 is the constant from Lemma 


example we consider the following three cases. We abbreviate 


6.3 


For this 


4 1124 

,( “ ):=,( “-3 ) = Ha!i7T 6 

Case 1 (unconstrained problem) In this case we set 


CL 


3 a 24 . 


lib OO, 

yb = -y a = oo. 


The values of Hp^H^a, r](a) and J(uh) for different values of a with choice A1 
for yo are given in Table [ 8 ] The findings are illustrated graphically in Figure [9] 
We see that Uh is a global minimum for all values of a since ||ph||z,s is less than 
its corresponding 77 (a). On the other hand, with choice A2 for y$ we can claim 
that Uh is a global minimum only for approximately a greater than 1 as it can 
be seen from Figure [lOj The numerical values are provided in Table [9] 

Case 2 (constrained control) In this case we consider constraints only on the 
control, we set 


u a = -5, 

u b = 5, 

Vb = -~Va = oo- 


Table 10 shows the values of ||p/i||l 6 , ?y(a) and J(uh) computed for different 
values of a with choice A1 for y 0 . The graphical illustration of these findings are 
shown in Figure |11| We see that Uh is a global minimum for a approximately 
greater than The numerical results associated with the choice A2 are 

given in Table |TT| and illustrated in Figure 12 In this case Uh is a global 
minimum for a approximately greater than 1 . 


Case 3 (constrained state) In this case we consider constrains only on the state, 
we set 


Ub — u a — oo, 

Va = - 1 , 
yb = 1 - 

The numerical findings associated with choice A1 are provided in Table [12] and 
illustrated in Figure [13] We see that Uh is a global minimum for all values of a. 
For the choice A2, the results are given in Table [13| and illustrated in Figure [14] 
We see that Uh is a global minimum only for a approximately greater than 1. 
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Table 7 Example [l] Case 4 : The values of ||p/i||i4, 77(a) and J(uh) for different 
values of a. 


a 

IbhIU 4 

77 (a) 

J i^h) 

1.0e-06 

1.961933031441e-04 

6.776197632762e-03 

2.143984056211e-01 

1.0e-05 

7.663887131231e-04 

1.606889689070e-02 

2.410556714493e-01 

1.0e-04 

2.844056064106e-03 

3.810535956559e-02 

2.890783107664e-01 

1.0e-03 

1.055630139945e-02 

9.036204771862e-02 

3.690000948128e-01 

1.0e-02 

2.397197977885e-02 

2.142821839497e-01 

4.449373232494e-01 

1.0e-01 

4.706175447556e-02 

5.081431366100e-01 

4.917394785652e-01 

1.0e+00 

4.818113594926e-02 

1.204997272869e+00 

4.991551130306e-01 

1.0e+01 

4.829535470702e-02 

2.857498848277e+00 

4.999153188201e-01 

1.0e+02 

4.830679945384e-02 

6.776197632762e+00 

4.999915299530e-01 

1.0e+03 

4.830794415727e-02 

1.606889689070e+01 

4.999991529760e-01 


Table 8 Example [ 2 ] Case 1 with choice A1 for yo : The values of ||ph||i6, 77(a) and 
J(u h ) for different values of a. 


a 

IKIk 6 

77 (a) 


1.0e-06 

1.179795342411e-04 

8.697974773247e-04 

3.663839269975e-03 

1.0e-05 

1.040717291260e-03 

2.498914960443e-03 

3.314332555914e-02 

1.0e-04 

6.486412414763e-03 

7.179344781194e-03 

1.967178952607e-01 

1.0e-03 

1.467650352720e-02 

2.062614866979e-02 

4.320253853445e-01 

1 . 0 e -02 

1.672495487678e-02 

5.925861229879e-02 

4.922543706340e-01 

1 . 0 e -01 

1.696149575588e-02 

1.702490943800e-01 

4.992144828609e-01 

1 . 0 e +00 

1.698552077353e-02 

4.891230660460e-01 

4.999213370331e-01 

1 . 0 e +01 

1.698792705311e-02 

1.405243150394e+00 

4.999921325890e-01 

1 . 0 e +02 

1.698816771892e-02 

4.037242258255e+00 

4.999992132478e-01 

1.0e+03 

1.698819178587e-02 

1.159893577654e+01 

4.999999213247e-01 


Table 9 Example [ 2 ] Case 1 with choice A2 for yo'. The values of ||p/i||x,6, 77(a) and 
J(v.h) for different values of a. 


a 

\\Ph\W 

77 (a) 

J(Uh) 

1.0e-06 

5.510426875132e-03 

8.697974773247e-04 

1.185192313978e+02 

1.0e-05 

1.587525748968e-02 

2.498914960443e-03 

1.331807740335e+02 

1.0e-04 

4.474831409415e-02 

7.179344781194e-03 

1.473322027953e+02 

1.0e-03 

1.039480114464e-01 

2.062614866979e-02 

1.584387338104e+02 

1 . 0 e -02 

2.428391864045e-01 

5.925861229879e-02 

1.626178840362e+02 

1 . 0 e -01 

3.493646725426e-01 

1.702490943800e-01 

1.642025836782e+02 

1 . 0 e +00 

3.554038724369e-01 

4.891230660460e-01 

1.644198119684e+02 

1 . 0 e +01 

3.560155910725e-01 

1.405243150394e+00 

1.64441976641 le+02 

1 . 0 e +02 

3.560769159456e-01 

4.037242258255e+00 

1.644441976184e+02 

1.0e+03 

3.560830499750e-01 

1.159893577654e+01 

1.644444197614e+02 
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( a ) \\Ph\\Le and 17(a) vs. a. 


(b) J(uh) vs. a. 



(c) The optimal state yh- 


(d) The adjoint state p^- 



(e) The optimal control Uh- 

Figure 9 Example [ 2 ] Case 1 with choice A1 for yo: The values of ||jht||_L6, 77 (a) and 
J{uh) vs. a. The optimal state yn, the optimal control Uh and the adjoint state ph 
for a = 10~ J . 
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(a) IIp/iIIlb and P(«) vs- a. 



(b) J(«h) vs. a. 



(c) The optimal state yh- 


(d) The adjoint state ph- 



(e) The optimal control uh- 

Figure 10 Example [ 2 ] Case 1 with choice A2 for yo: The values of ||ph||ne, ry(a) 
and J{u h ) vs. a. The optimal state yh, the optimal control un and the adjoint state 
Ph for a = 1. 
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Table 10 Example [ 2 ] Case 2 with choice A1 for yo: The values of ||ph|| L 6, 77(a) and 
J{uh) for different values of a. 


a 

IKIU 6 

77 (a) 

J iV'h) 

1.0e-06 

1.613825290585e-02 

8.697974773247e-04 

4.507855415302e-01 

1.0e-05 

1.613824266503e-02 

2.498914960443e-03 

4.508885528139e-01 

1.0e-04 

1.613816501602e-02 

7.179344781194e-03 

4.519051721159e-01 

1.0e-03 

1.615565078678e-02 

2.062614866979e-02 

4.612661359991e-01 

1 . 0 e -02 

1.672495487678e-02 

5.925861229879e-02 

4.922543706340e-01 

1 . 0 e -01 

1.696149575588e-02 

1.702490943800e-01 

4.992144828609e-01 

1 . 0 e +00 

1.698552077353e-02 

4.891230660460e-01 

4.999213370331e-01 

1 . 0 e +01 

1.698792705311e-02 

1.405243150394e+00 

4.999921325890e-01 

1 . 0 e +02 

1.698816771892e-02 

4.037242258255e+00 

4.999992132478e-01 

1.0e+03 

1.698819178587e-02 

1.159893577654e+01 

4.999999213247e-01 


Table 11 Example^ Case 2 with choice A2 for yo: The values of ||p?i||l6, 77(a) and 
J(u h ) for different values of a. 


a 

IKIU 6 

77 (a) 

J(Uh) 

1.0e-06 

3.456649663660e-01 

8.697974773247e-04 

1.636832040856e+02 

1.0e-05 

3.456649990198e-01 

2.498914960443e-03 

1.636833073745e+02 

1.0e-04 

3.456663172695e-01 

7.179344781194e-03 

1.636843379602e+02 

1.0e-03 

3.456602557101e-01 

2.062614866979e-02 

1.636944643396e+02 

1 . 0 e -02 

3.457537810584e-01 

5.925861229879e-02 

1.637855203878e+02 

1 . 0 e -01 

3.494672249476e-01 

1.702490943800e-01 

1.642029145907e+02 

1 . 0 e +00 

3.554038724369e-01 

4.891230660460e-01 

1.644198119684e+02 

1 . 0 e +01 

3.560155910725e-01 

1.405243150394e+00 

1.64441976641 le+02 

1 . 0 e +02 

3.560769159456e-01 

4.037242258255e+00 

1.644441976184e+02 

1.0e+03 

3.560830499750e-01 

1.159893577654e+01 

1.644444197614e+02 


Table 12 Example [ 2 ] Case 3 with choice A1 for yo: The values of ||Ph||i 6 , 77(a) and 
J{u h ) for different values of a. 


a 

IKIU 6 

77 (a) 

J(Uh ) 

1.0e-06 

1.293594798095e-04 

8.697974773247e-04 

6.247856764953e-02 

1.0e-05 

8.673961098825e-04 

2.498914960443e-03 

8.936458658379e-02 

1.0e-04 

5.421978025542e-03 

7.179344781194e-03 

2.033602173575e-01 

1.0e-03 

1.467650352720e-02 

2.062614866979e-02 

4.320253853445e-01 

1 . 0 e -02 

1.672495487678e-02 

5.925861229879e-02 

4.922543706340e-01 

1 . 0 e -01 

1.696149575588e-02 

1.702490943800e-01 

4.992144828609e-01 

1 . 0 e +00 

1.698552077353e-02 

4.891230660460e-01 

4.999213370331e-01 

1 . 0 e +01 

1.698792705311e-02 

1.405243150394e+00 

4.999921325890e-01 

1 . 0 e +02 

1.698816771892e-02 

4.037242258255e+00 

4.999992132478e-01 

1.0e+03 

1.698819178587e-02 

1.159893577654e+01 

4.999999213247e-01 
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(a) Uphllie and 77 (a) vs. a. 
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(b) J(u h ) vs. a. 





(c) The optimal state yh- 




(e) The optimal control uh- 
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(f) The control active sets (Uh = 5 inside 
the green circles and Uh = — 5 inside the 
blue ones). 


Figure 11 Example [ 2 ] Case 2 with choice A1 for yo : The values of ||ph||n 6 , 77 (a) 
and J{u h ) vs. a. The optimal state yh, the optimal control Uh, the control active 
sets, and the adjoint state ph for a = 10 ~ 3 . 
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( a ) \\Ph\\Le and 17 (a) vs. a. 


(b) J(uh) vs. a. 



(c) The optimal state yh- 


(d) The adjoint state ph- 



(e) The optimal control uh- 

Figure 12 Example [2] Case 2 with choice A2 for yo: The values of ||ph||i6, ri(a) 
and J{u h ) vs. a. The optimal state yh, the optimal control un and the adjoint state 
Ph for a = 1. 
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(c) The optimal state yn- 


(d) The adjoint state ph- 



A md A 



(e) The optimal control Uh- 


(f) The multipliers /rjj (in blue) and p b h 
(in green). 


Figure 13 Example [ 2 ] Case 3 with choice A1 for yo : The values of ||ph||n6, p(a) 
and J(uh) vs. a. The optimal state ph, the multipliers p^,p b h , the optimal control uh 
and the adjoint state ph for a = 10~ 5 . 
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( a ) IIp/i IIlb and 17(a) vs. a. ^ J ( Uh ) vs ' a 



(c) The optimal state yh- 


(d) The adjoint state ph- 



(e) The optimal control uh- 

Figure 14 Example [2] Case 3 with choice A2 for yo: The values of ||ph||ne, »?(a) 
and J{u h ) vs. a. The optimal state yh, the optimal control un and the adjoint state 
p h for a = 1. 
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6 Appendix 


Lemma 6.1 We have for a, b > 0, A, y > 0 that 

a x b» < AA ^ (a + b) x+ ». 

~ (A + y) x +» K 

Proof: Apply Young’s inequality xy < -p x p + ^ y® x, y > 0, -p + ^ = 1 

to P = ^ \ ^ , Q = ^ ^ and a: = (Pa) p , y = (Qb) Q . 

A y 

Lemma 6.2 Suppose that Assumption 1 holds. Then we have for a,i£l 
^ <j>'(ta + (1 — t)b) — <j)'(b) dt\ < \a — b\L r ( f ft (ta + (1 — t)b) dt) r , 


where 


r — 1 , r=± 


=- 


dt 


Proof: We start by noticing that 

f 4>'(ta + (1 — t)b) — ()>'(b) dt = f f (f"(ra + (1 — r)b)(a — b) dr 

Jo Jo Jo 

= (a — b) f (1 — t)(j>" (ta + (1 — t)b) dt. 
Jo 

Therefore, taking the absolute value and using Assumption 1 we get 

r 1 r 1 i 

(j>'(ta + (1 — tfb) — 4>'(b) dt <\a — b\M / (1 — t)(j>'(ta + (1 — t)b) r dt 

o Jo 

fi i 

< \a-b\M\\l-t\\ L r'( 0A) ( / <j>'(ta + (1 - t)b) dtY, 


where A + p = 1. It is easy to see that 

Denoting M\\l — t\\ Lr >( 0 l) by A r completes the proof. ■ 


Theorem 6.3 (Gagliardo—Nirenberg interpolation inequality^ 

2 

For 2 < q < oo we define 9 = 1 - as well as 

q 


GN q \= sup 

/e-F/HR 2 ),/^) 


II /|| Lq(M 2 ) 

II/IIl 2 (r 2 ) liv/ll 


e 

L 2 ( R 2 ) 


Then GN q < C q := min(G^, Cf \ Cf >), where 
C<p = (9C 2ae )-\ ifq> 4; 


( 6 . 1 ) 
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C (2) = - 1 - 

9 V0 e (i-0) 


= (2nB(l, 2 ^ — -)) 9/2 k B (—^y, (6.2) 

i_ e V v ) 2 q )> By 2 + 29 J y ’ 


-\ _ 2 OO 

c? = <-) % n 


2i 


;=2 23 +q - 2 


2^+2 -q 
I 2 3 q 


(6.3) 


Here, 


C 2 , s = 2 1 / s (^) (s 1)/s (27ri3(?,3-^)) 1/2 , l<s<2; C' 2>1 =2 Vtt; 


B (M) = M>0 

r(a + o) 

, / N ( P \1/P(P' \-l/p' 1 . 1 1 

k B{P) = (x~) (z-) , - + - = 1. 

v 27r v 27r 7 p p 


Proof: The bounds (6.1) and (6.2) can be found in the paper m by Veling. 

(1.7)]. The estimate 

is QS (1.42),(1.43)], 


We remark that GN q = A 2 where A 2 0 is defined in m 


(6.1) is [T2! (1.31)] (note that 9 > - O q > 4), while (6.2 


where the latter bound has been proved by Nasibov in m 


Let us now turn to the proof of (6.31. To begin, we claim that for all k £ No 

k+l 21 ai+a-t i k - 


3 —2 


^|/ll^l|V/|| i2 % (6.4) 


where 


q k = 2~ k (q + 2(2 k -1)). 


The inequality clearly holds for k = 0. Suppose that (6.41 is true for some 
k £ No- We infer from Theorem 1 in [5] for the case d = 2 that 

ll/IU 2 *’ < ^ll/llip+i l|V/||?, 2 , 1 < p < oo. (6.5) 

Here, 

/i = ( i , ( p r 1)2 )h^)*( „ r(i, >.. 'i f wM, o = 2(p r 1) . y = p+1 


47T 


2 y ’ V r(p-1)' 


4p 


P~ 1 


Using the formula for y and observing that T(y) = (y — l)T(y — 1), the expression 
for A can be simplified to 

1 e 71 -U I e i 

A =(-r(^r~ 2p - 


We apply (6.51 for p = -q k and obtain 


where 


\l« < Aii/rr hv/ii^, 

L2 9 fe t 


+ „ Qk -2 


( 6 . 6 ) 


A =(-)>( and 9 = 


^Qk 
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and 9 = 1 


Qk +1 


Since -Qk + 1 = Qfe+i we find that 




Qk 


which, inserted into ( 6.61 yields 


i < (IyC-^y -n/n^ +1 iiv/r. 


(6.7) 


Using the induction hypothesis we infer 




i \ gfc 


fe+i 

n< 

i=2 


2-? 


2^ + g - 2 


2 J '+2-fj l-5iL+(95i 

I 2j? ll/llU+x 9 IIV/II^ 5 '■ 


Elementary calculations show that 


1 0 Ok 
2 ' q ' 2 q 




(i — 0)— 


^ /i _ Qk+l \ 

2 '■ g 
2 k +2 


' 2 fc + 2 + g — 2 ; 
9fc+i 


2 fc + 2 _[_ 2 — g 
| 2 fc + 2 g 


= 1 - 


9fc+i 

q 


which implies (6.4) for k + 1. The result now follows by sending k —> oo in (6.41 
and by observing that lim g*, = 2. 

/c —yoo 


35 

















References 


[1] N. Arada, E. Casas, and F. Troltzsch. Error estimates for a semilinear elliptic 
control problem. Computational Optimization and Applications 23:201-229 
( 2002 ) 

[2] E. Casas. Error estimates for the numerical approximation of semilinear el¬ 
liptic control problems with finitely many state constraints. ESAIM Control 
Optimisation and Calculus of Variations 8:345-374 (2002). 

[3] E. Casas and F. Troltzsch. Second order optimality conditions and their role 
in PDE control. Jahresbericht der Deutschen Mathematiker-Vereinigung 
DOI 10.1365/sl3291-014-0109-3 (2014). 

[4] K. Deckelnick and M. Hinze. A finite element approximation to elliptic 
control problems in the presence of control and state constraints. Hamburger 
Beitrage zur Angewandten Mathematik 2007-01 (2007). 

[5] M. Del Pino, J. Dolbeault. Best constants for Gagliardo-Nirenberg inequali¬ 
ties and applications to nonlinear diffusions. J. Math. Pures Appl. 81:847-875 
( 2002 ). 

[6] M. Hintermiiller, K. Ito, and K. Kunisch. The primal-dual active set strategy 
as a semismooth Newton method. SIAM Journal on Optimization, 13:865- 
888 ( 2002 ). 

[7] M. Hinze. A variational discretization concept in control constrained op¬ 
timization: The linear-quadratic case. Computational Optimization and 
Applications, 30:45-61 (2005). 

[8] M. Hinze and A. Rosch. Discretization of optimal control problems. Inter¬ 
national Series of Numerical Mathematics 160:391-431 (2011). 

[9] M. Hinze, R. Pinnau, M. Ulbrich, S. Ulbrich. Optimization with pde con¬ 
straints. Mathematical Modelling: Theory and Applications, Volume 23, 
Springer (2008). 

[10] S.M. Nasibov On optimal constants in some Sobolev inequalities and their 
application to a nonlinear Schrodinger equation. Soviet. Math. Dokl. 40:110- 
115 (1990), translation of Dokl. Akad. Nauk SSSR 307:538-542 (1989). 

[11] I. Neitzel, J. Pfefferer, and A. Rosch. Finite element discretization of state- 
constrained elliptic optimal control problems with semilinear state equation. 
SIAM Journal on Control and Optimization, to appear (2015). 

[12] M. Ulbrich. Semismooth Newton methods for operator equations in func¬ 
tion spaces. SIAM Journal on Optimization, 13(3):805-841, 2002. 

[13] E.J.M. Veling. Lower bounds for the infimum of the spectrum of the 
schrodinger operator in and the sobolev inequalities. JIPAM. Journal of 
Inequalities in Pure & Applied Mathematics [electronic only], 3, Art. 63 
( 2002 ). 


36 



Table 13 Example [ 2 ] Case 3 with choice A2 for yo: The values of ||ph||i, 6 , 77 (a) and 
J(u h ) for different values of a. 


a 

IKIU 8 

?y(a) 

J(Uh) 

1.0e-06 

1.139290773221e-03 

8.697974773247e-04 

1.525635040951e+02 

1.0e-05 

8.200728224157e-03 

2.498914960443e-03 

1.536016384574e+02 

1.0e-04 

2.474482888749e-02 

7.179344781194e-03 

1.559116076253e+02 

1.0e-03 

9.716506658549e-02 

2.062614866979e-02 

1.600204462920e+02 

1 .0e-02 

1.800129125912e-01 

5.925861229879e-02 

1.627566303073e+02 

1 .0e-01 

3.493646725426e-01 

1.702490943800e-01 

1.642025836782e+02 

1 .0e+00 

3.554038724369e-01 

4.891230660460e-01 

1.644198119684e+02 

1 .0e+01 

3.560155910725e-01 

1.405243150394e+00 

1.64441976641 le+02 

1 .0e+02 

3.560769159456e-01 

4.037242258255e+00 

1.644441976184e+02 

1.0e+03 

3.560830499750e-01 

1.159893577654e+01 

1.644444197614e+02 
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